Pairing Symmetry and Upward Curvature of Upper Critical Field in Superconducting 
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The origin for the upward curvature of the upper critical field (H C 2) observed in hydrate cobaltate 
Nao.35Co02 ■ yHbO is investigated based on the microscopic gap equation. It is shown that the 
observed upward curvature results from the transition between two different pairing symmetries 
that occur on different energy bands. Furthermore, different pairing symmetries involved in the 
transition results in different upward curvatures. By considering transitions among all lowest possible 
pairing symmetries, it is found that the transition of the pairing symmetry from s-wave at low fields 
to d x 2_ y 2 + id xy at high fields is the best fit to the experimental data. Our results provide an 
important clue to the understanding of the superconductivity in hydrate cobaltate. 



PACS numbers: 73.20.-r, 73.20.At, 73.21.Hb 

Since the discovery of superconductivity in hydrate 
cobaltate Nao.ssCoCVyEbOp]], extensive theoretical and 
experimental studies have been devoted to elucidate the 
mechanism of superconductivity. To unravel the mecha- 
nism, identifying the underlying pairing symmetry would 
be the first step. The two-dimensional triangle lattice 
formed by C0O2 provides an alternative lattice symmetry 
to the square CuC>2 lattice in high Tc materials and has 
led to many proposals for unconventional pairing symme- 
tries. For singlet pairing, the lowest possible unconven- 
tional symmetry without breaking rotational symmetry 
is d x 2_ y 2 ± id xy , and mechanism based on the correla- 
tion effects for such symmetry has been proposed Im- 
possible spin-triplet /-wave pairing was also proposed 
based on ferromagnetic fluctuations 0, 0, 0]. Exper- 
imentally, however, data reported show contradictory 
conclusions 8] , indicating the fragility of the supercon- 
ductivity in this system. Furthermore, there are evi- 
dences indicating that two pairing symmetries may be 
involved in this system 0, [To, EH • The issue of whether 
the pairing symmetry is conventional or not is still un- 
settled and needs to be clarified. 

In this paper, we shall focus on the data of the upper 
critical fields, measured by the specific heat. The specific- 
heat technique probes the bulk properties of the samples 
and has been proved to be a powerful tool for inves- 
tigating the pairing state of many superconductors |l2j| . 
Specifically, for the hydrate cobaltate Nao.asCoCVyH^O, 
upward curvatu re ( a kink structure in the slope) of H C 2 
was observed @, [nj • Similar structure was also observed 
in early studies of high Tc materials .13]. Based on the 
Ginzburg-Landau theory, Joynt[3] attributed the up- 
ward curvature to the transition between two different 
pairing symmetries with different critical temperatures. 
However, in a later investigation based on microscopic 
formulation of the gap equation [3], negative results 
were found, indicating that the upward curvature is not 



due to mixing of two pairing order parameters. The rea- 
son why two approaches give different results lies in the 
fact that in the Ginzburg-Landau theory, H C 2 oc T — T c 
and phenomenologically, both T c and slopes of H& are of- 
ten chosen arbitrarily. If larger slope of H c2 is chosen for 
smaller T c , H C 2 of two pairing symmetries near their T c 's 
essentially shows the intersection of two straight lines. A 
kink in the slope thus arises and the upward curvature 
can be easily simulated. In real materials, however, the 
slope of H C 2 and T c both depend on microscopic details 
and are not independent from each other. In fact, in the 
Gorkov's microscopic derivation of the Ginzburg-Landau 
equation [nij], the slope is proportional to m*T c /ep with 
to* being the effective mass of the electron and ep be- 
ing the Fermi energy. For a single band, to* ftp are the 
same for different pairing symmetries, hence smaller T c 
goes with smaller slope, in the opposite trend adopted 
in the Ginzburg-Landau equation. Hence in this case, 
joining two pairing symmetries with different T c would 
not yield the upward curvature. This picture essentially 
explains why the upward curvature is not reproduced in 
the calculation of Kim et al. (HI for high Tc materi- 
als. Then, what is the origin for the upward curvature in 
hydrate colaltate ? It is known that multi-orbitals near 
the Fermi surface might be involved for the occurrence 
of superconductivity [1 7| . Hence to* / 'ep can no longer be 
treated as a fixed parameter for different pairing sym- 
metries if different pairings occur on different bands. In- 
deed, our calculation in below shows that upward curva- 
ture feature can result from the two-band calculation in 
which different values of m* jtp are assumed for different 
energy bands where different pairing symmetry occurs. 
Furthermore, mixing of different pairing symmetries re- 
sults in coupling of the ground state to different Landau 
levels in the presence of magnetic fields and causes differ- 
ent upward curvatures. By direct comparison of experi- 
mental data with H C 2 obtained by mixing of the lowest 
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possible pairing symmetries, it is possible to pin down 
the pairing symmetries of the system. We find that the 
upward curvature observed in the experimental data is 
mostly consistent with a transition of the pairing sym- 
metry from s to d x 2_ y 2 + id xy . 

We start by first considering pairing that occurs in two 
bands. The two-dimensional BCS-like Hamiltonian can 
be written as 

2 

H = E E ^c^dk^ 

+ E E Vkik'jCktfC-kilC-k'jlCk'jh (!) 

which was first investigated by Suhl, Matthias, and 
Walker [18]. Here ki and k2 are wave vectors on two 
Fermi surface sheets indexed by 1 and 2. The electron- 
electron interaction Vkik'j and H^k;, are the intrasheet 
contributions and V^ki, is an intersheet contribution. In 
general, the superconducting pairing symmetry is differ- 
ent for different band. Hence we can write the projection 
of the interaction in the form Vkik' = ^^(ki)^* (k^), 
V k2k > = V^(k 2 )<^(ky, and V klk > = ViiaQLjfaQ^). 

Here a and /3 are indices for the pairing symmetry. (f> a (k) 
is an operator that projects out the corresponding pair- 
ing symmetry a. For the lowest possible pairings, we 
shall consider possible mixing of p and / waves for triplet 
pairing and mixing of s and d for singlet pairing. For the 
triangle lattice, assuming that the rotational symmetry 
is not broken, the appropriate p and d pairing amplitudes 
are p x ± ip y and d x 2_ y 2 ± id xy . The corresponding pro- 
jection operators are </> s (k) = 1, p ±i p (k) = - k x ± ik y , 
4>d±id{k) = k 2 x -k 2 y ± 2ik x k y , and </>/(k) = k% - 3k x k%. 
Note that there are two possible /-waves; however, since 
they are related by a rotation of tt/6 and it suffices to 
consider one of them (see below for details). Following 
Kim et al. 15], the real-space linearized gap equation in 
thepresence of the magnetic field V x A can be written 

as nil 

A a (R) = V a J2 I dT^{<j>)^ a {(t>)k x {v,uj)A a {K) 
+VjJ2 I drfc(<p)fo(<j>)K 2 (r,u)Ap(R), (2) 

and 

+ViJ2 [ dxfa {<t>)M4>)Ki (r,w)A (R). (3) 

Here R represents the position for the center of mass of 
the Cooper pair, r is the displacement of the Cooper pair 
and <j) is the corresponding angle of r. A's are pairing 



amplitudes in real space and <j) a {<l>) are the projection 
operators in real space: <fr s {<f>) — 1) < f > p±ip( < t>) = e± ^ 
<t ) d±id{4>) = e ±2l( ^, and </>/(0) = cos 3 ^ — 3 cos 4> sin 2 
cos3(f>. Ki(r,u}) is the kernel operator, given by 



Ki(r, w) = Kf(r, w) exp r • (y R 
with 

K9(r,Lj) = k B TN i (0y 
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where kpi, VFi and -/Vj(O) are the Fermi wavenumber, 
Fermi velocity and the two dimensional density of state 
for the ?th band. In the absence of Vj, Eqs.© and © 
decouples, and their solutions for constant A and A = 
yield relations between T c and V a 19] 
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where i = a or (3 with T % c being the correspond- 
ing transition temperature in zero field and r, — 
^\4>i(4>)\ 2 - Dividing Eqs.@ and © by V a and 
Vp respectively and using Eq.©, V a and Vp can be 
eliminated. Following Ref.[l9j, if we further adopt di- 
mensionless variables, ti = T/TJ., r = pj y / 2eH/hc and 
hi = 2eH /hc{hvFi/2iikBTl) 2 , the gap equations become 
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(7) 
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where = Vj/Vi and 7 = Np/N a . The operators K are 
given by 
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where a™ = 2\L^~ VFm and a = V + ^eA) x ± 



2^ 

i(V + 2ie J 4) B ]. 

Eqs.([7]) and ([8]) are the governing equations for the 
situation when two pairing symmetries occur on different 
energy bands. However, it also covers the case when the 
two pairing symmetries occur in the same energy band. 
In that case, one simply sets 7 = 1, £j = 1, and drop the 
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index i. Eqs.([7]) and |(5J) can be solved by expanding A Q 
and A.p in terms of the set of eigenfunctions \ipn) °f the 
two-dimensional Schrodinger equation in the presence of 
A = (0,Hx,0) 



and 
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n=0 



(10) 



where i — a or i — (3. Note that the eigenfunctions 
\ip n ) are essentially the Landau levels. The operators 
K couple different Landau levels. The detailed cou- 
pling is determined by the bi-product projection opera- 
tor, (j)^((j))(j) m ((/>), in Eq.©, which, after being integrated, 
selects the correct combinations of a + and a~ that sur- 
vive. The selected combination of a + and a - then de- 
termines how A % n couple. For the mixing of s-wave and 
d + id, {A%} couples with {A^^}; while for the mix- 
ing of s-wave and d — id, {A!^} couples with {A„t^ f}. 
On the other hand, for the mixing of p ± ip and /-wave, 
because 4>f(<f) contains both e 3 ^ and e _3l< ^, in addition 
to the coupling between A J ^ tlp and A^, there are also 
couplings among Since the bi-product projection 

operator in K/f, <j)*f{(j))4>f{(j)), contains e ±t6 ^ , A^ cou- 
ples with A^ +6 for each n; while because the bi-product 
projection operators in off-diagonal K a p (a ^ j3), con- 
tains e ±z4 * and e ±l2 ^, we found that {A^} couple with 
{A p n XJ} and {A*} couple with {A£$}. The coupling 
coefficients are most conveniently expressed in terms of 
the following functions [l5| 
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where i is the index for the pairing symmetry, L n {x) 
is the Laguerre polynomial and L^ k is the associated 
Laguerre polynomial in Rodrigues representation. For 
triplet pairing, wc find that recursion relations for the 
mixing of /-wave with p x ± ip y are given by 
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(6,n-6K_ 6 + ln-+/?n A{ 
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-^7\/^i?« P±lp (T4, n ± 4)^±f = 0, (13) 
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On the other hand, for singlet pairing, recursion relations 
for the mixing of s-wave with d x 2_ y 2 ± id xy are given by 
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( T 2,n±2)A d n ± ± f = 0, 
(15) 



\n- + f3i)A d n ±ld + £ >d±ldl 



a s (±2,n T 2)A i nT2 =Q. 

(16) 



When applying the above analysis to the calculation 
of H C 2 with mixing of two pairing symmetries, a and 
(3, one assumes that the pairing symmetries are associ- 
taed with different T c 's with T c a > if. In other word, 
the pairing symmetry a is the stable bulk pairing state 
at low fields while j3 is the stable pairing state at high 
fields. Therefore, one starts from A a = Aq\iPq). Mix- 
ing to the other symmetry, (3, then couples Aq to A@ 
with n > 1. Since {A s n } couples with {A d ^ 2 d }, this 
analysis implies that one can only have the transitions 
from high-temperature s-wave to low-temperature d + id 
or from high-temperature d — id to low-temperature s- 
wave. On the other hand, for the mixing of p ± ip 
and /-wave, because </>/(</>) contains both e 3 ^ and e~ 3 ^, 
both transitions from high-temperature p ± ip to low- 
temperature /-wave and from high-temperature /-wave 
to low-temperature p±ip are possible. Note that there 
are two possible / — Avaves, k x — <^k x k y and k y — 3k y k x . 

Since the two /-waves are related by exchanging k x and 
k y which simply exchanges p x +ip y andpz— ip y , it suffices 
to consider one of them. 

When solving H C 2, it is important to note that there 
are many eigenvalues H(T) satisfying Eqs.® and ([5]) 
and only the largest one defines H c2 . To compare the 
calculated H c2 with experimental data, one needs to fix 
scales of temperatures and magnetic fields. The tran- 
sition temperature T c of the most stable bulk pairing 
state determines the temperature scale. On the other 
hand, the scale of magnetic fields can be fixed by the 
data points with lower magnetic fields. The remaining 
parameters are the ratio of Fermi velocities, the ratio of 
T c , £j and 7. At this point, it is important to note that 
according to Eqs. p3|) - (fl6|) . for different mixing scheme, 
different Landau levels are mixed. As a result, differ- 
ent mixing scheme results in different upward curvature. 
To find the best fit to the experimental data, after fix- 
ing scales of temperatures and magnetic fields, we vary 
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the remaining parameters to find the best. We find that 
the transition from s-wave at low fields to d + id at high 
fields is the best fit to the data. In Fig. 1, we show nu- 
merical results of H c2 for the transition from s-wave at 
low fields to d + id at high fields in comparison with ex- 
perimental data obtained by specific- heat measurement. 
The fitting parameters are the ratio of Fermi velocities 
Vp/vp-" = 0.8, T?+ id /T° = 0.5, & = 4.3, U+id = 0.89, 
and 7 = 1.1. These values are in reasonable regime. The 
close fitting ro the experimental data clearly shows that 
singlet pairing dominates in hydrate cobaltat, which is 
also consistent with recent NMR data[20]. Furthermore, 
it implies that two energy bands are involved and sup- 
ports results based on LDA calculations [l?} where two 
bands constructed from the three Co tid orbitals inter- 
sect the Fermi level. 




3.0 4.0 
T(K) 



FIG. 1: The comparison of experimental data of -ff C 2 0] with 
numerical results based on the transition from s-wave at low 
fields to d + id at high fields. 

In conclusion, we have investigated the origin for the 
upward curvature of the upper critical field {H C 2) ob- 
served in hydrate cobaltate Nao.3sCo02 • J/H2O. Analysis 
based on the microscopic gap equation shows that the 
observed upward curvature results from the transition 
between two different pairing symmetries that occur on 
different energy bands. Furthermore, it is found that the 
transition of the pairing symmetry from s-wave at low 
fields to d x 2_ y 2 + id xy at high fields is the best fit to the 
experimental data. 
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